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1.  Introduction 

In  this  paper  we  develop  lower  and  upper  bounds  for 
arrowhead  matrices.  A  matrix  Q  G  IRmxm  is  called  an 
arrowhead  matrix  if  it  has  a  form  as  follows: 

D  c~ 


where  D  G  a  diagonal  matrix,  c  is  a  vector 

in  IR.”1^1,  and  b  is  a  real  number.  Here  the  superscript 
“t”  signifies  the  transpose.  The  arrowhead  matrix  Q  is 
obtained  by  bordering  the  diagonal  matrix  D  by  the  vector 
c  and  the  real  number  b.  Hence,  sometimes  the  matrix 
Q  in  (1)  is  also  called  a  symmetric  bordered  diagonal 
matrix.  In  physics,  arrowhead  matrices  have  been  used  to 
describe  radiationless  transitions  in  isolated  molecules  [1] 
and  oscillators  vibrationally  coupled  with  a  Fermi  liquid  [2]. 
Numerically  efficient  algorithms  for  computing  eigenvalues 
and  eigenvectors  of  arrowhead  matrices  were  discussed  in 
[3].  The  properties  of  eigenvectors  of  arrowhead  matrices 
were  studied  in  [4],  and  as  an  application  of  their  results,  an 
alternative  proof  of  Cauchy’s  interlacing  theorem  was  given 
there.  The  existence  of  arrowhead  matrices  was  investigated 


recently  in  [5-8]  such  that  the  constructed  arrowhead 
matrix  has  the  pregiven  eigenvalues  and  other  additional 
requirements. 

Our  motivation  to  study  lower  and  upper  bounds  of 
arrowhead  matrices  is  from  Kung  and  Suter’s  recent  work  on 
the  hub  matrix  theory  [9]  and  its  applications  to  multiple- 
input  and  multiple  output  (MIMO)  wireless  communication 
systems.  A  matrix,  say  A,  is  a  hub  matrix  with  m  columns  if  its 
first  m—  1  columns  (called  nonhub  columns)  are  orthogonal 
to  each  other  with  respect  to  the  Euclidean  inner  product 
and  its  last  column  (called  hub  column)  has  a  Euclidean 
norm  greater  than  any  other  columns.  Subsequently,  it  was 
shown  that  the  Gram  matrix  of  A,  that  is,  Q  =  A* A,  is 
an  arrowhead  matrix  and  its  eigenvalues  could  be  bounded 
by  the  norms  of  the  columns  of  A.  As  pointed  out  in 
[9-11],  the  eigenstructure  of  Q  determines  the  properties 
of  wireless  communication  systems.  This  motivates  us  to 
reexamines  these  bounds  of  the  eigenvalues  of  Q  and  makes 
them  sharper.  In  [9],  the  hub  matrix  theory  is  also  applied 
to  the  MIMO  beamforming  problem  by  comparing  k  of  m 
transmitting  antennas  with  the  largest  signal-to-noise  ratio, 
including  the  special  case  where  k  =  1  which  corresponds 
to  a  transmitting  hub.  The  relative  performance  of  resulting 
system  can  be  expressed  as  the  ratio  of  the  largest  eigenvalue 
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of  the  truncated  Q  matrix  to  the  largest  eigenvalue  of  the 
Q  matrix.  Again,  it  was  previously  shown  that  these  ratios 
could  be  bounded  by  the  ratios  of  norms  of  columns  of  the 
associated  hub  matrix.  Sharper  bounds  will  be  presented  in 
Section  4. 

The  well-known  result  on  the  eigenvalues  of  arrowhead 
matrices  is  the  Cauchy  interlacing  theorem  for  Hermitian 
matrices  [12].  We  assume  that  the  diagonal  elements  dj, 
j  =  1, 2, . . . ,  m  —  1,  of  the  diagonal  matrix  D  in  ( 1)  satisfy 
the  relation  d\  <  d2  <  ■  ■  ■  <  dm-i-  Let  Ai,  A2, . . . ,  Am  be  the 
eigenvalues  of  Q  arranged  in  increasing  order.  The  Cauchy 
interlacing  theorem  says  that 

^1  —  d\  <  A2  ^  —  ■  ■  ■  —  dm— 2  —  Am_i  <  dm-\  <  Am.  (2) 


Lemma  1.  Let  Q  E  K.mxm  be  an  arrowhead  matrix  of  the  form 
(1),  where  D  =  diag(di,d2»--->^m-i)  £  ||(m-i)x(m-i)j  b  E  R, 
andc  =  (ci,C2,...jCm-i)  E  [Rm_1.  Then 

m—  1  m- 1  2m—  1 

det(Xl-Q)  =  (X-b)Y\(X-dk)-  X \cj[Y\(X-dkl 

k= 1  ;'  =  1  k=l 

kfj 

(6) 

The  proof  of  this  result  can  be  found  in  [5,  13]  and 
therefore  is  omitted  here. 

When  the  diagonal  matrix  D  in  ( 1 )  is  a  zero  matrix,  the 
following  result  is  followed  from  Lemma  1 . 


When  the  vector  c  and  the  real  number  b  in  (1)  are  taken 
into  consideration,  a  lower  bound  of  Aj  and  an  upper  bound 
of  Am  were  developed  by  using  the  well-known  Gershgorin 
theorem  (see,  e.g.,  [3,  12]),  that  is, 


m—  1 


Am  <  max|iii  +  |ci  dm-i  +  \cm-i\,b+  X  kilj,  (3) 

Ai  >  mint  d\  -  |ci  |,. . . ,  dm-i  -  |cm_i|,fc-  X  Id  \-  (4) 


i=  1 


Accurate  bounds  of  eigenvalues  of  arrowhead  matrices 
are  of  great  interest  in  applications  as  mentioned  before. 
The  main  results  of  this  paper  are  presented  in  Theorems 
1 1  and  12  for  the  upper  and  lower  bounds  of  the  arrowhead 
matrices.  It  is  also  shown  in  Corollary  13  that  the  resulting 
bounds  are  tighter  than  in  (2),  (3),  and  (4). 

The  rest  of  the  paper  is  outlined  as  follows.  In  Section  2, 
we  will  introduce  notation  and  present  several  useful  results 
on  the  eigenvalues  of  arrowhead  matrices.  We  give  our 
main  results  in  Section  3.  In  Section  4,  we  revisit  the  lower 
and  upper  bounds  of  the  ratio  of  eigenvalues  of  arrowhead 
matrices  associated  with  hub  matrices  and  wireless  com¬ 
munication  systems  [9],  and  subsequently,  we  make  those 
bounds  shaper  by  using  the  results  in  Section  3.  In  Section  5, 
we  compute  the  bounds  of  arrowhead  matrices  using  the 
developed  theorems  via  three  examples.  Conclusions  are 
given  in  Section  6. 


Corollary  2.  Let  Q  E  M.mxm  be  an  arrowhead  matrix  having 
the  following  form: 


Q 


(7) 


where  c  is  a  vector  in  IRm  1  and  b  is  a  real  number.  Then  the 
eigenvalues  of  Q  are 


b  -  ^b2  +  4||c||2  b  +  -\jb2  +  4||c||2 

2  >  —  2  > 

A ;(Q)  =  0,  for  i  =  2, . . . ,  m  -  1. 

(8) 


Proof.  By  using  Lemma  1,  we  have 

det(A7  -  Q)  =  Am~2(A2  -bX-  ||c||2).  (9) 

Clearly,  A  =  0  is  a  zero  of  det(A/  -  Q)  with  multiplicity  m  —  2. 
The  zeros  of  the  quadratic  polynomial  A2  -  bX  -  ||c||2  are 
( b  -  fb2  +  4||c||2)/2  and  ( b  +  fb2  +  4||c||2)/2,  respectively. 
This  completes  the  proof.  D 


In  what  follows,  a  matrix  Q  having  a  form  in  (7)  is  called 
a  special  arrowhead  matrix.  The  following  corollary  (also,  see 
[3])  is  a  direct  result  from  Lemma  1. 


2.  Notation  and  Basic  Results 

The  identity  matrix  is  denoted  by  I.  The  notation 
diag(fli,  fl2>  •  •  •  >  (in)  represents  a  diagonal  matrix  whose  diag¬ 
onal  elements  are  a\ ,  <22, . . . ,  an.  The  determinant  of  a  matrix 
A  is  denoted  by  det(A).  The  eigenvalues  of  a  symmetric 
matrix  A  E  RnXn  are  always  ordered  such  that 

Ai  (A)  <  A 2(A)  <  ■  ■  ■  <  A „{A).  (5) 

For  a  vector  a  E  K”,  its  Euclidean  norm  is  defined  to  be 
Ha||  :=  I  a,- 12. 

The  first  result  is  about  the  determinant  of  an  arrowhead 
matrix  and  is  stated  as  follows. 


Corollary  3.  Let  Q  be  an  m  X  m  arrowhead  matrix  given  by 
(1),  where  D  =  diag(di,d2,---.dm-i)  E  [rC1”— i)x(m— 1)^  b  E  K, 
and  c  =  (ci,c2,...,  cm- 1)  E  IRm_I.  Let  us  denote  the  repetition 
of  the  number  dj  in  the  sequence  {d,-}”  j  by  kj.  If  kj  >  2,  then 
dj  is  the  eigenvalue  of  Q  with  multiplicity  kj  —  1. 

Proof.  When  the  integer  kj  >  2,  the  result  follows  from 
Lemma  1  since  (A  —  dj)ki  l  is  a  factor  of  the  polynomial 
det(A7  -  Q).  Id 

Corollary  4.  Let  Q  be  an  m  x  m  arrowhead  matrix  given 
by  { 1),  where  D  =  diag(di,  d2,. . .  ,dm-i)  £  |j^(m-i)x(m-i)) 
b  E  M,  and  c  =  (ci,  C2, . . . ,  cm- 1)  e  r^1.  Suppose  that  the 
last  k  >  2  diagonal  elements  dm_^,  dm-/t+ 1,  ■  ■  ■  >  dm- 1  of  D  are 


2 


EURASIP  Journal  on  Advances  in  Signal  Processing 


3 


identical  and  distinct  from  the  first  m-k  -  l  diagonal  elements 
d\,di,...,  d„,-i t_i  ofD.  Define  a  new  matrix 


[di 


Cl 


Q:= 


Cl 


d-m-k-l  Cm—k-1 

d-tn—k  Cm-k 
Cm—k— 1  Cm-k  b 


(10) 


with  Cj  =  Cj  for  j  =  1,2, . . .  ,m  —  k  —  1  and  cm-k  = 
■\jllf=m-k  \cj\2-  Then  the  eigenvalues  ofQ  are  that  ofQ  together 
with  dm-k  with  multiplicity  k  -  1. 


entries  of  the  vector  c  in  ( 1 )  are  nonzero.  The  reason  for  this 
assumption  is  the  following.  Suppose  that  C;,  the  jth  entry  of 
c,  is  nonzero,  it  can  be  easily  seen  from  Lemma  1  that  A  -  dj 
is  a  factor  of  det(A/  -  Q);  that  is,  dj  is  one  of  eigenvalues  of 
Q.  The  remaining  eigenvalues  of  Q  are  the  same  as  those  of 
a  matrix  which  is  obtained  by  simply  deleting  the  jth  row 
and  column  of  Q.  In  summary,  for  any  arrowhead  matrix, 
we  can  find  eigenvalues  corresponding  to  repeated  values  in 
D  or  associated  with  zero  elements  in  c  by  inspection. 

In  this  paper,  we  call  a  matrix  Q  in  (1)  irreducible  if 
the  diagonal  elements  di,  d2, . . . ,  dm-i  of  Q  are  distinct  and 
all  elements  of  c  are  nonzero.  By  using  Corollary  4  and  the 
above  discussion,  this  arrowhead  matrix  can  be  reduced  to 
an  irreducible  one. 


Proof.  Since  numbers  dm-ib  dm-/t+i,. . .  ,  d„,-i  are  identical 
and  distinct  from  numbers  di,  d2, . . .  ,dm-k- 1>  we  have 


m-1 


n  (^  -  di) 

i=i 


/  \ 

m-k 

n(A- di) 
1=1 

Vti  ) 


(A  -dm.k)k~\ 


j  < m-k -  1 , 


Remark  5.  In  [4,  9],  Hermitian  arrowhead  matrices  are 
considered;  that  is,  it  allows  that  c  in  the  matrix  Q  of  the  form 
(1)  is  a  vector  in  We  can  directly  construct  many  (real 
symmetric)  arrowhead  matrices  denoted  by  Q  from  Q.  The 
diagonal  elements  of  these  symmetric  arrowhead  matrices 
are  the  exactly  same  as  those  of  Q.  The  vector  c  in  Q  could 
be  chosen  as 


m-1  2m 

I  | Cj  I  UiX-di) 

j=m—k  i=l 

i*j 

(m-1  2\  / m-k—1  \ 

x  | Cj |  n  (A-d,o \(x-dm-k)k-\ 

j=m—k  J  \  i=  1  / 

(ID 


c=  (±|Ci|,±|C2l,...,±|cm-i|).  (13) 

In  such  a  way,  there  are  2m~'  such  symmetric  arrowhead 
matrices.  Because  det(A7  -  Q)  =  det(Af  -  Q)  by  Lemma  1, 
every  such  symmetric  arrowhead  matrix  Q  has  the  identical 
eigenvalues  with  Q.  This  is  the  reason  why  we  just  consider 
the  eigenvalues  of  real  arrowhead  matrices  in  this  paper. 


By  (6)  in  Lemma  1,  we  have 
det(A7  -  Q) 


The  following  well-known  result  by  Weyl  on  eigenvalues 
of  a  sum  of  two  symmetric  matrices  is  used  in  the  proof  of 
our  main  theorem. 


m-1  m-1  2m_^ 

=  (\-b)Y\tt-di)-  X  \  ci  |  m-di) 

i=  1  ;  =  1  i=  1 

it) 


Theorem  6  (Weyl).  Let  F  and  G  be  two  m  x  m  symmetric 
matrices.  Let  us  assume  that  the  eigenvalues  of  F,  G,  and  F  +  G 
have  been  arranged  in  increasing  order.  Then 


( 

m-k 

(x-b)Yl(x-di)- 

i=i 

V 


m-k 

1 


;=i 


.m-k-l 

n 


\ 

(A -dd 


2=1 

/ 


(A -dm-k)k~l 


=  det(A7  -  Q)  ■  (A  -dm-k)k~\ 


(12) 


Clearly,  if  A  is  an  eigenvalue  of  Q,  then  A  is  either  an 
eigenvalue  of  Q  or  dm-k-  Conversely,  dm-k  is  an  eigenvalue 
of  Q  with  multiplicity  k  -  1  and  the  eigenvalues  of  Q  are  that 
of  Q.  This  completes  the  proof.  EUl 

By  using  Corollaries  3  and  4,  to  study  the  eigenvalues  of 
Q,  we  may  assume  that  the  diagonal  elements  d\,  d2,  ■  ■  ■ ,  dm-i 
of  Q  are  distinct  when  we  study  the  eigenvalues  of  Q  in 
( 1 ) .  Since  eigenvalues  of  square  matrices  are  invariant  under 
similarity  transformations,  we  can  without  loss  of  generality 
arrange  the  diagonal  elements  to  be  ordered  so  that  d\  < 
da  <  ■  ■  ■  <  dm- 1.  Furthermore,  we  may  assume  that  all 


Xj(F  +  G)  <  A,(F)  +  Aj_i+m(G), 

for  i>  j. 

(14) 

A j(F  +  G)  >  Xj(F)  +  A;-j+i(G), 

for  i  <  j. 

(15) 

Proof.  See  [14,  page  62]  or  [12,  page  184] 

□ 

To  apply  Theorem  6  for  estimating  eigenvalues  of  an 
irreducible  arrowhead  matrix  Q,  we  need  to  decompose  Q 
into  a  sum  of  two  symmetric  matrices  whose  eigenvalues  are 
relatively  easy  to  be  computed.  Motivated  by  the  structure 
of  the  arrowhead  matrix  and  the  eigenstructure  of  a  special 
arrowhead  matrix  (see,  Corollary  2),  we  write  Q  into  a  sum 
of  a  diagonal  matrix  and  a  special  arrowhead  matrix. 

To  be  more  precisely,  let  Q  e  [Rmxm  be  an  irreducible 
arrowhead  matrix  as  follows: 


Q 


D 


(16) 


3 


4 

where  dm  e  RL,  D  =  diag(di,  d2, . . . ,  dm_i)  with  0  <  d\  < 
d2  <  ■  ■  ■  <  dm_ 1  <  dm,  and  c  is  a  vector  in  [Rm_1.  For  a  given 
p  e  [0, 1],  we  write 

Q  =  £  +  S,  (17) 


Therefore,  we  can  use  Theorem  6  to  give  estimates  of  the 
eigenvalues  of  Q  via  those  of  E  and  S.  To  number  the 
eigenvalues  of  E,  we  introduce  the  following  definition. 

Definition  7.  For  a  number  p  e  [0, 1],  we  define  an  operator 
Tp  that  maps  a  sequence  {d; }J=l  satisfying  0  <  d\  <  d2  < 

■  ■  ■  <  dm_ 1  <  dm  to  a  new  sequence  {dj}f=l  :=  Tp({di}"‘=1) 
according  to  the  following  rules:  if  pdm  <  d\,  then  di  :=  pdm 
and  dj+i  :=  dj  for  j  =  1  ,...,m  —  1;  if  pdm  >  dm_i,  then 
dj  :=  dj  for  j  =  —  1  and  dm  :=  pdm;  otherwise, 

there  exists  an  integer  j0  such  that  dj0  <  pdm  <  dj0+ 1,  then 

dj  .  dj  for  j  1,. . . ,  j0,  dj0+\  .  pdm ,  and  dj+i  .  dj  for 

=  j0  +  -  1. 

Theorem  8.  Let  Q  e  IRmxm  be  an  irreducible  arrow¬ 
head  matrix  having  a  form  of  (16),  where  D  = 
diag(di,d2>- ■  ■  >dm_i)  with  0  <  d\  <  d2  <  ■  ■  ■  <  dm-i  <  dm, 
and  c  is  a  vector  in  IRm_1.  For  a  given  p  E  [0,1],  define 
{di}njl=l  :=  rp({d,};m=1).  Then,  one  has 


min 

di 

+  t,  d2,  dm  +  sj, 

if  j  =  1. 

Aj(Q)  <  - 

min 

dj 

+ 1,  dj+i  j, 

if  2  <  j  <  m 

-  1, 

dm  + 

t, 

if  j  =  m, 

(19) 

d\  + 

s, 

ifi  = 

1, 

O 

IV 

max 

[dj 

-l,  dj  +  sj , 

if  2< 

j  * 

m  -  1, 

max 

d\ 

+  t ,  dm_  i ,  dm  I-  s 

>  fj  = 

m, 

(1  -  p)dm  -  ^(l-p)2d2m+4\\c\\2 

s  =  - - - , 

2 

(21) 

=  (l-p)dm+V(l-p)2^+4||c||2 

2 

Proof.  For  a  given  number  p  e  [0, 1],  we  split  the  matrix  Q 
into  a  sum  of  a  diagonal  matrix  E  and  a  special  arrowhead 
matrix  S  according  to  ( 1 7) ,  where  E  and  S  are  defined  by  ( 1 8 ) . 
Clearly,  we  know  that 

A  j(E)  =  dj  (22) 
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for  j  =  1, 2, . . . ,  m.  By  Corollary  2,  we  have 

Ai(S)=s,  Am(S)  =  t,  Aj(S)  =  0,  for  j  =  2,...,m  -  1, 

(23) 

where  s  and  t  are  given  by  (21). 

Upper  Bounds.  By  ( 14)  in  Theorem  6,  we  have 
Xj(Q)  <  Xi(E)  +  X 

m+j-  i(S )  (24) 

for  all  i  >  j.  Clearly,  for  a  given  j, 

X ;(Q)  <  minjAiCE)  +A  m+j—i  (S)}-  (25) 

More  precisely,  since  {d/}™  j  is  monotonically  increasing,  s  < 
0,  and  t  >  0,  we  have 

Ai(Q)  <  minjdi  +  t,  d2, . . . ,  dm  +  sj 

=  min  jdi  +  t,  d2,  dm  +  s| , 

A ;(Q)  <  min | dj  +  f,dj+i,...,dm}  =  minjdj  +  t,dj+ 1| 

(26) 

for  j  =  2, . . . ,  m  -  1,  and 

A m(Q)  ^  A m(E)  +  X m(S)  =  dm  +  t.  (27) 

In  conclusion,  (19)  holds. 

Lower  Bounds.  By  ( 15)  in  Theorem  6,  we  have,  for  a  given  j, 
Aj(Q)  >  max|A,(E)  +  Aj_;+i(S)|.  (28) 

Hence, 

Ai(Q)  >  Ai(£)  +  Ai(S)  =  d\  +s, 

Aj(Q)  >  maxjdj  +  s,dj~ i,. . . , di  |  =  maxjdj  +  s, dj_i| 

(29) 

for  j  =  2, . . . ,  m  —  1,  and 

Am(Q)  >  max jdm  +  s,  dm-i,. . . ,  d2,  d\  +  t| 

=  max|dm  +  s,  d„,_i,  d\  +  f  j. 

As  we  can  see  from  Theorem  8,  the  lower  and  upper 
bounds  of  the  eigenvalues  for  Q  are  functions  of  p  e  [0, 1  ]  for 
the  given  irreducible  matrix  Q.  In  other  words,  the  bounds 
of  eigenvalues  vary  with  the  number  p.  Particularly,  when  we 
choose  p  being  the  ending  points,  that  is,  p  =  0  and  p  =  1, 
we  can  give  an  alternative  proof  of  interlacing  eigenvalues 
theorem  for  arrowhead  matrices  (see,  e.g.,  [12,  page  186]). 
This  theorem  is  stated  as  follows. 
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Theorem  9  (Interlacing  eigenvalues  theorem).  Let  Q  e  RmXm 
be  an  irreducible  arrowhead  matrix  having  a  form  in  (16), 
where  D  =  diag(di,  d2,. . .  ,dm-\)  with  0  <  d\  <  da  <  ■  ■  ■  < 
dm-  i  <  dm,  and  c  is  a  vector  in  Let  the  eigenvalues  of  Q 
be  denoted  by  withX\  <  X2  <  ■  ■  ■  <  Xm.  Then 

Ai  <  <4  <  A2  —  d2  ^  ^  dm—2  —  Am_j  <  dm-\  <  Am. 

(31) 


Prop/!  In  Theorem  8,  the  upper  bounds  of  the  eigenvalues  of 
Q  in  ( 19)  are  determined  by  dj,  j  =  1,2,...  ,m,  and  s  and  t 
in  (21).  They  can  be  viewed  as  functions  of  p  in  [0, 1],  That 
is,  the  upper  bounds  of  the  eigenvalues  of  Q  are  functions  of 
p  in  the  interval  [0,1].  Therefore,  we  are  able  to  find  optimal 
bounds  of  the  eigenvalues  of  Q  by  choosing  proper  p.  The 
upper  bounds  on  XjiQ)  for  j  =  1,  2  <  j  <  m  -  1,  and  j  =  m 
in  (34)  are  discussed  separately. 


Proof.  By  using  ( 19)  with  p  =  0  in  Theorem  8,  we  have 
A  j  <  dj  for  j  =  1,2,. . . ,  m  -  1.  By  using  (20)  with  p  =  1 
in  Theorem  8,  we  obtain  A  j  >  dj- 1  for  j  =  2,3  ,...,m. 
Combining  these  two  parts  together  yields  our  result.  □ 

The  proof  of  the  above  result  shows  that  we  could 
have  improved  lower  and  upper  bounds  for  each  eigenvalue 
of  an  irreducible  arrowhead  matrix  by  finding  an  optimal 
parameter  p  in  [0,1].  Our  main  results  will  be  given  in  the 
next  section. 

3.  Main  Results 

Associated  with  the  arrowhead  matrix  Q  in  Theorem  8,  we 
define  four  functions  f,  i  =  1, 2, 3, 4,  on  the  interval  [0, 1]  as 
follow: 

flip)  :=\((l-p)dm-^(l-p)2d]n+4\\c\\2), 

flip)  :=  \((l-p)dm  +  ^l-p)2dl+m\2),  (32) 

hip)  :=  pdm  + flip), 
hip)  ■■=  pdm+  flip)- 

Obviously, 

s  =  fiip),  t  =  f2ip),  (33) 


Upper  Bound  ofXiiQ).  From  ( 19),  we  have 

Ai (S)  <  min jdi  +  f,d2,(im +  sj,  (35) 

where  <4,  s,  and  t  are  functions  of  p  on  the  interval 
[0,1].  In  this  case,  we  consider  p  in  the  following  four 
subintervals.  [0,  d\! dy^\,  [ d ] / ,  di/ d J ,  [ di/ d^^ ,  d^  . \ l d,,, ] , 
and  [dm-i/dm,  1],  respectively.  For  p  e  [0,di/dm],  we  have 
d]  +  t  =  hip),  di  =  di,  and  dm  +  s  =  flip).  For  p  e 
[d]/dm, di!dm),  we  have  d\  +  t  =  d\  +  flip),  di  =  pdm,  and 
dm+  s  =  dm- 1  +  flip).  For  p  e  [d2/ dm,dm-\l dm],  we  have 
di+t=di+  flip),  d2  =  d2,  and  dm  +  s  =  dm- 1  +  flip).  For 
p  e  [dm-i/dm,  1],  we  have  d\  +  t  =  di  +  flip),  di  =  d2,  and 
dm  +  s  =  hip).  Hence 


min  d2  =  di, 

pe[  0,1] 


where  s  and  t  are  given  by  (21). 

The  following  observation  about  monotonicity  of  func¬ 
tions  f,  i  =  1, 2, 3, 4,  is  simple,  but  quite  useful  as  we  will  see 
in  the  proof  of  our  main  results. 

Lemma  10.  The  functions  /]  and  f2  both  are  decreasing  while 
h  and  h  ar£  increasing  on  the  interval  [0,1]. 

The  proof  of  this  lemma  is  omitted. 

Theorem  11.  Let  Q  be  an  irreducible  arrowhead  matrix 
defined  by  (16)  and  satisfying  all  assumptions  in  Theorem  8. 
Then  the  eigenvalues  of  Q  are  bounded  above  by 


XjiQ)  < 


minjdi.dm-!  +  /i(^J‘)}> 
dj, 


dm- 1  + 


if  j  =  1, 

if  2  <  j  <  m  —  1, 


min(di  +  t 

pev  ' 


fiiO),  ifV  = 

/  /L  \ 

if  V  = 

dm—  1 


dl+f2it )’ 


0 

d\  d2 
d, 


rfl  +  /2(^f)’  ifV“ 

di  +  fiil),  if  V  = 


dy] 

'dm  r 

dy\ 

'dm  J’ 

dm- 1  "] 

>  —j  > 

Clm  J 


(36) 


Since  0  >  fi(di/dm)  >  fiid2/dm)  >  fi{dm-i/dm),  /2(0)  > 
dm  >  du  and  f2id2/dm)  >  fiidm-i/dm)  >  /2(  1)  >  0,  we  have 

Ai(Q)  <  mm^dudm-i  (37) 


if  j  =  m. 


(34) 


5 


Upper  Bound  of  XjiQ),  for  2  <  j  <  m- 1.  From  (19),  we  have 
Xj(Q)  <  minjdj  +  t,dj+ 1|.  (38) 
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In  this  case,  we  consider  p  lying  in  the  following  four  subin¬ 
tervals.  [0,  dj— \/ djyf  ] ,  [ dj — ] / dm,  ci j / dm  | ,  [dj/dm ,  dj+i/dm ] ,  and 
[dj+}/dm,  1],  respectively.  For  p  e  [0,dj-i/dm],  we  have  dj  + 
t  =  dj-\  +  f2{p)  and  dj+i  =  dj.  For  p  e  [dj-i/dm,dj/dm],  we 
have  dj  +  t  =  /4(p)  and  dj+ 1  =  dj.  For  p  e  [dj/dm,  dj+\/dm], 
we  have  dj  +  t  =  dj  +  flip)  and  d,+1  =  pdm.  For  p  e 

[dj+i/dm,  1],  we  have  dj  +  t  =  dj  +  flip)  and  dj+i  =  dj+\. 
Hence 


min  dj+i  =  dj, 
pe[o,i] 


min(di  +  t)  = 


dj- i  +  /2 
dj-1  +  /2 
dj  +  /2 


dj—i 
dm 
dj- 1 


2j+l 


+  /2(1), 


if  V  = 


if  V  = 


ifV  = 


ifV  = 


0, 


dj-i 


dj- 1  dj 

dm  dm 

dj  dj+i 
dm  dm 


dj- 1 


,1 


(39) 


Therefore, 


jdj-i  +  /2^ 


dj- 1 


Aj(Q)  <  min-j  dj- 1  +  /2  (  -y—  ] ,  d;-  [. 

Since  d;_i  +  /2(dj_r 
we  get 

Aj(Q)  <  d;-. 

Upper  Bound  of\m{Q).  From  (19)  we  have 
Am  (  Q )  —  ^m  3*  t . 


(41) 


(42) 


For  p  e  [0,  d„,_i/dm],  we  have  dm  +  t  =  dm~i  +  flip)  while 
for  p  e  [dm_i/dm,  1],  we  have  dm  +  t  =  /4(p): 


min  (dm  +  f)  = 


peV 


Hence, 


dm-!  +/2r' 
V  dm 

j  ,  r  l  dm  - 1 


II 

> 

i  i 

© 

3  i 

i _ i 

II 

i  i 

a.  f»* 

3  7 

I—1 

-T"  l _ l 

dm- 1  \ 

dm  / 

(' 

This  completes  the  proof. 


□ 


Theorem  12.  Let  Q  be  an  irreducible  arrowhead  matrix 
defined  by  (16)  and  satisfying  all  assumptions  in  Theorem  8. 
Then  the  eigenvalues  of  Q  are  bounded  below  by 


A  AQ)> 


max-jdj-i,  dj  +  fi{^~ 

dl+f2{t)’ 


if  j  =  1, 
if  2  <  j  <  m  — 
if  j  =  m. 

(45) 


Proof.  In  Theorem  8,  the  lower  bounds  of  the  eigenvalues 
of  Q  in  (20)  are  determined  by  dj,  j  =  1,2 ,...,m,  and  s 
and  t  in  (21).  As  we  did  in  Theorem  12,  the  lower  bounds 
of  the  eigenvalues  of  Q  are  functions  of  p  in  the  interval 
[0, 1].  Therefore,  we  are  able  to  find  optimal  bounds  of  the 
eigenvalues  of  Q  by  choosing  proper  p.  The  discussion  is 
given  for  j  =  1,  2  <  j  <  m  -  1,  and  j  =  m  in  (45),  separately. 


Lower  Bound  o/Ai(Q).  From  (20),  we  have 
Ai  (Q)  >  d\  +  s. 


(46) 


In  this  case,  we  consider  p  lying  in  the  following  two 
subintervals:  [0, d\/dm\  and  [d\/dm,  1].  For  p  e  [0, di/dm], 
di  +s  =  flip).  For  p  e  [dj/dm,  1],  we  have  di  -l-s  =  d\  +  flip). 
Hence 


(47) 


(40) 

maxfdi  +  s)  =  - 

di  +  fi 

'H 

Him  ' 

av  = 

k  j-l 

Um  J 

—  dm  >  dj , 

peV  1  > 

di  +  fi 

K  Him  ' 

av  = 

ViA 

-  HLm  J 

It  leads  to 


Lower  Bound  ofXiiQ).  From  (20),  we  have 
A2(Q)  >  max|di,d2 +s|. 


(48) 


(49) 


In  this  case,  we  consider  p  lying  in  the  following  three 
subintervals:  [0,di/dm],  [di/dm,d2/dm],  and  [d2/dm,l].  For 
p  e  [0,di/dm],  we  have  dj  =  pdm,  d2  +  s  =  di  +  flip).  For 
p  e  [di/dm,d2/dm],  we  have  di  =  d\,  d2  +  s  =  flip).  For 
p  e  [di/dm,  1],  we  have  d\  =  d2  and  d2  +  s  =  d2  +  flip). 
Hence, 


max(d2  +  s)  = 

peV  '  ' 


maxdi  =  di, 

pev 


di  +  fdO), 
di 


Hi)-  ■ 

Hi)’  ‘ 


d2  +  fi  (  ~r~ 
di  +  fi  (  -f- 


HV  = 

if  V  = 

HV  = 


0,dl\ 

d\  d-2  ~| 
dm  dm  J 

f.il. 

um  J 


(50) 
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These  lead  to 

A2(Q)  >  maxjdi,d2  +  fi  |-  (51) 

Lower  Bound  ofXj(Q),  3  <  j  <  m  -  1.  From  (20),  we  have 
Aj(Q)  >  max|dj_i, dj  +  s|.  (52) 

In  this  case,  we  consider  p  lying  in  the  following  three  subin¬ 
tervals.  [0,  dj— 2/dm [ ,  [ dj — 2/ dm,  dj  —  1  / dm\y  [ dj  ] / dm ,  dj/ dm  J , 
and  [dj/dm,  1].  For  p  G  [0, dj-2/dm],  we  have  dj_i  =  dj_2 
and  dj  3-  s  =  dj- 1  +  /i(p).  For  p  G  [dj-2/dm,  dj-\/dm], 
we  have  dj_i  =  pdm  and  dj  +  s  =  dj-\  +  fi{p).  For  p  G 
[dj_i/dm,dj/dm],wehavedj_i  =  dj- 1  and  dj+s  =  flip).  For 

p  G  [dj/dm,  1],  we  have  dj- 1  =  dj_i  and  dj  3-  s  =  dj  3-  /i(p). 
Hence 


maxd,_i  =  dj-i, 


max  (dj  3-  5)  = 


dj- 1  3-  MO), 

if  y  = 

o,a-p 

dm 

j 

di— 2 

di— 1 

if  y  = 

j  a 

dm 

j _ 

dm 

(  dj\ 

\  um  / 

r  , 

di  1 

if  y  = 

j  j 

dm  dm  _ 

(  dj  \ 

+  /'.j  ’ 

\  dm  / 

r  dj  1 

if  y  = 

dm 

(53) 


Since  dj- 1  >  dj_i  +  /i(0)  >  dj_i  +  /i(dj_2/d„,),  we  have 


Aj(Q)  >  max]dj„i,dj  +  fi  I  '-j- 

V  dn 


(54) 


Lower  Bound  of\m(Q).  From  (20),  we  have 

Am(Q)  >  maxjdi  +  t,dm-\,dm  3-sj-.  (55) 

In  this  case,  we  consider  p  lying  in  the  following  three  subin¬ 
tervals.  [0,di/dm],  [ d| / dm , dm  2/ dm ] ,  [dm—2/dm, dm—i/dm\, 
and  [d„,_i/dm,l].  For  p  G  [0,di/dm],  we  have  d\  +  t  = 
f^ip),  dm- 1  dm—2>  dm  3"  s  dm-\  3-  f\ip).  For  p  G 

[ di/ dm,  dm— 2/ dm ]  ,  we  have  di  3- 1  di  -t  (p),  dm— 1  dm— 2, 
dm  3-  s  —  dm— 1  3~  ftip)-  For  p  G  [dm— 2/ dm , dm—  i/ dm ] ,  we  have 
di  3-  f  =  di  3-  /2(p),  dm-i  =  pdm,  dm  +  s  =  dm- 1  3-  /i(p).  For 


p  £  [dm— 1/ dm,  1],  we  have  di  3- 1  —  d\  3-  flip ),  dm-i  —  dm-i, 
dm  3-  s  =  flip).  Hence 


max  dm — 1  —  dm— 1, 

pHO.l] 


max(dm  3-  s) 

p£V  v  > 


max(di  3-  f) 
pev v  ’ 


dm-i  +  MO), 


dm-i+  fi  (  -j-  )> 


7  .  r  (  dm- 2 

“m- 1+  /l  I  “j - 

\  wm 

dm  +  /l(l), 

dl+f2it )’ 

di+/2(dl 
di  +■  /a  ^ 


d 

dm-2 


d\  3- 


r  ( dm— l 

\  am 


ifV  = 

if  V  = 

,  ifV  = 

ifV  = 

if  y  = 

if  y  = 

if  y  = 
if  y  = 


dm  J 

d\  dm—2 
dm  dm 

dm- 2 


\ 


dm- 1  1 

um  J 


—  ll 

7  >  1  > 

ww  J 

o,Al 

dm  J 

^1  dm— 2 


dm  dm 

dm- 2  dm 
dm  dy, 


f]' 


-4 

1  J 


(56) 


Since  0  >  /i(0)  >  /i(di/dm)  >  /i(dm-2/dm)  and/2(di/dm)  > 
fi(dm  J dm )  >  /2(dm-i/dm),  we  have 

Am(Q)  >  max|dm-i,dm  3- /i(l),di  3- /2(^-)  j.  (57) 


Since  di  3-  /2(di/dm)  =  Md\/dm)  >  MO)  >  dm,  we  get 

Am(Q)>d1+/2(^).  (58) 

This  completes  the  proof.  □ 

Corollary  13.  Let  Q  be  an  irreducible  arrowhead  matrix 
defined  by  (16)  and  satisfying  all  assumption  in  Theorem  8. 
Then  upper  and  lower  bounds  of  the  eigenvalues  of  Q  obtained 
by  Theorems  11  and  12  are  tighter  than  those  given  by  (2),  (3), 
and  (4). 


Proof.  Since 


min|di,dm-i  3-/i^^-d^|  <  di,  (59) 


then  the  upper  bound  for  the  eigenvalue  Ai(Q)  given  by  (34) 
in  Theorem  11  is  tighter  than  that  by  (2).  The  upper  bounds 
for  the  eigenvalues  A jiQ),  j  =  2, . . . ,  m  -  1 ,  provided  by  (34) 
in  Theorem  1 1  are  the  same  as  those  by  (2). 
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Note  that  0  <  d\  <  ■  ■  ■  <  dm-\  <  dm;  the  right-hand  side 
of  (3)  with  b  =  dm  is 


I  d\  +  | Ci  | ,  dm~ i  +  \cm~i  I,  b  +  X  kil 


m- 1 

=  dm  +  I  | . 
2=1 


(60) 


Since  ||c||  <  YULi  kil,  dm  +  ||c|l  =  Ml),  and  dm-X  + 
fi{dm-i/dm)  =  Mdm-i/dm),  we  have 


m- 1 

dm  I  Cj  | 
2=1 


(61) 


and  then  the  upper  bound  of  A„,(Q)  from  (34)  in  Theorem  11 
is  tighter  than  that  from  (3). 

Now  we  turn  to  the  lower  bounds  of  Aj(Q).  Since 


(62) 


4.  Hub  Matrices 

Using  the  improved  upper  and  lower  bounds  for  the 
arrowhead  matrix,  we  will  now  examine  their  applications  to 
hub  matrices  and  MIMO  wireless  communication  systems. 
The  concept  of  the  hub  matrix  was  introduced  in  the  context 
of  wireless  communications  by  Kung  and  Suter  in  [9]  and  it 
is  reexamined  here. 

Definition  17.  A  matrix  A  £  K"xm  is  called  a  hub  matrix,  if  its 
first  m—  1  columns  (called  nonhub  columns)  are  orthogonal 
to  each  other  with  respect  to  the  Euclidean  inner  product 
and  its  last  column  (called  hub  column)  has  its  Euclidean 
norm  greater  than  or  equal  to  that  of  any  other  columns.  We 
assume  that  all  columns  of  A  are  nonzeros  vectors. 

We  denote  the  columns  of  a  hub  matrix  A  by 
«i,  fl2,  ■  ■  ■ ,  flm-  Vectors  a\,  ai,.  ■ .  ,  am-i  are  orthogonal  to  each 
other.  We  further  assume  that  0  <  ||flil|  <  || a2 II  <  ■  ■  ■  < 
||am||.  In  such  case,  we  call  A  an  ordered  hub  matrix.  Our 
interest  is  to  study  the  eigenvalues  of  Q  =  A‘A,  the  Gram 
matrix  A.  In  the  context  of  wireless  communication  systems, 
Q  is  also  called  the  system  matrix.  The  matrix  Q  has  a  form 
as  follows: 


for  j  =  2, . . . ,  m  —  1  and 


d\  + 


^  dm  dyi 


(63) 


llfltll2 


\\a2\\2 


Q  = 


{d\ ,  CLm  ) 
(#2,  tt-m  ) 


we  know  that  the  lower  bounds  for  the  eigenvalues  A j(Q), 
j  =  2, provided  by  (45)  in  Theorem  12  are  tighter 
than  those  by  (2).  Ip 


II  dm—  1 II 

{fl-mi  d\)  {dm>  CI2)  ’  "  *  {d  m , Am— 1 ) 


1 ,  ) 

II  Q-m  II 

(66) 


Remark  14.  When  c  in  (16)  is  a  zero  vector,  by  using 
Theorems  11  and  12,  we  have  dj  <  A (Q)  <  dj,  that  is, 
A(Q)  =  dj.  In  this  sense,  the  lower  and  upper  bounds  given 
in  Theorems  1 1  and  12  are  sharp. 

Remark  15.  When  Q  in  Theorems  11  and  12  has  size  of 
2x2,  the  upper  and  lower  bounds  of  its  each  eigenvalue  are 
identical.  Actually,  from  Theorems  11  and  12  we  have 


Clearly,  Q  is  an  arrowhead  matrix  associated  with  A. 

An  important  way  to  characterize  properties  of  Q  is  in 
terms  of  ratios  of  its  successive  eigenvalues.  To  this  end,  the 
ratios  are  called  eigengap  of  Q  which  are  defined  [9]  to  be 


EGi(Q) 


Am-(i-l)  (Q) 
Am-i(Q) 


(67) 


d_ 1  +  /1  -  MQ)  ^  minjdi,di  +  /1  (  .  )  j> 

(64) 

dl+/2(f)  ^A2(Q)<di+/2(|). 

Clearly,  we  have 

Ai(Q)  =  di+/i(|),  A2(Q)  =  d1+/2(|).  (65) 

This  can  be  verified  by  calculating  the  eigenvalues  Q  directly. 


for  i  =  1,2, ...,m  —  1.  Following  the  definition  in  [9],  we 
define  the  ith  hub-gap  of  A  as  follows: 

HG,(A)  =  (68) 

I  am-i  II 

for  i  =  1,2, . . . ,  m  —  1. 

The  hub-gaps  of  A  will  allow  us  to  predict  the  eigen- 
structure  of  Q.  It  was  shown  in  [9]  that  the  lower  and  upper 
bounds  of  EGi  (Q)  [9]  are  given  by  the  following: 


Remark  16.  For  the  lower  bound  of  the  smallest  eigenvalue  HGi(A)  <  EGi(Q)  <  (HGi(A)  +  1)HG2(A).  (69) 

of  an  arrowhead  matrix,  no  conclusion  can  be  made  for  the 
tightness  of  the  bounds  by  using  (4)  and  (45)  in  Theorem  12. 

An  example  will  be  given  later  (see  Example  22  in  Section  5). 
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These  bounds  only  involve  nonhub  columns  having  the  two 
largest  Euclidean  norms  and  the  hub  column  of  A.  Using 
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the  results  in  Theorems  11  and  12,  we  obtain  the  following 
bounds: 


/4(||ai||2/||fl„,||2) 


I  ttm-1 I 


<  EGi(Q)  < 


M 

( II  ^m- 1 

llVlkJI2) 

maxj 

|  II &m— 

2 II 2 >  /3I 

[\\amM\2/\\am\\2)} 

(70) 


Obviously,  these  bounds  are  not  only  related  to  two  nonhub 
columns  with  the  largest  Euclidean  norms  and  the  hub 
column  of  A  but  also  related  to  the  nonhub  column  having 
the  smallest  Euclidean  norm  and  interrelationship  between 
all  nonhub  columns  and  the  hub  column  of  A.  As  we 
expected,  the  lower  and  upper  bounds  of  EGi(Q)  in  (70) 
should  be  tighter  than  those  in  (69).  To  prove  this  statement, 
we  give  the  following  lemma  first. 

Lemma  18.  Let  cq,  a2, . . . ,  am  be  the  columns  of  a  hub  matrix 
A  with  0  <  ||fli||  <  ||a2ll  <  ■  ■  ■  <  ||flm-ill  ^  ll«mll-  Then 


Mp)  >  llflmll2  for  p  e  (0,1], 


<  llflmll2  +  lla„,-il|2 


(71) 

(72) 


Proof.  From  Lemma  10,  we  know,  for  p  E  (0, 1], 

,,v  lkml|2+Vll«J|4+4||c||2  ..  ||2  (73) 

Mp)  >  Mo)  = - ^ - >  llflmll  ,  [/i> 


where  ||c||2  =  YT=i  I  (a,-,  am)  |2.  The  inequality  (71)  holds.  By 
the  definition  of  /4,  showing  the  inequality  (72)  is  equivalent 
to  proving 


|c||2  <  ||am||2||am-1||2. 


(74) 


This  is  true  because 


Proof  We  first  need  to  show 

II  am  II 2  f^Uif/UmW2) 

Clearly,  this  is  true  because  of  (71).  Next  we  need  to  show 


(76) 


M 

( II  &m—  1 

ll2/||am||2] 

1 

max] 

|  II &m— 

2 II  2>  /3I 

(ik.-illVilaJI2)} 

(77) 


+  1 


O-m- 1  I 


I  fl-m-  ill  /II  O-m—2  I 
To  this  end,  it  is  suffice  to  prove 


p-1  <  llflmll2  +  \\am-i\'2 
dm 


This  is  exactly  (72).  The  proof  is  complete. 


(78) 

n 


The  lower  bound  in  (70)  can  be  rewritten  in  terms  of  the 
hubgap  of  A  as  follows: 


/4(||flil|2/||a,„||2)  i 

— - =  ;HG  r(Q) 


I  O-m  - 1  I 


1  + 


1  + 


4 II  ell 


”\  II  O-n 


■  (79) 


The  upper  bound  in  (70)  can  be  rewritten  in  terms  of  the 
hubgap  of  A  as  follows: 


M 

^  II  Mm- 1 

Il2/Ilamll2) 

1 

max] 

|  II &m— 

2 12,  /3I 

(  km-lll2/  |flm  I2)} 

/4(||flm-1||2/||am||2) 


I  ^m-2  I 


=  -(HGi(A)  +  1)HG2(A) 


llaJI2  > 


m- 1 


s 


2 


> 


m- 1 


s 


(75) 


The  first  inequality  of  above  is  from  the  orthogonality  of  a;-, 
j  =  1, . . . ,  m  -  1  while  the  second  inequality  is  from  ||fli  ||  < 

1 1 a2 1 1  <  ■  ■  ■  <  ||am_i  ||.  This  completes  the  proof.  jjjj 


The  following  result  holds. 

Proposition  19.  Let  Q  in  (66)  be  the  arrowhead  matrix 
associated  with  a  hub  matrix  A.  Assume  that  0  <  || a4 1|  < 
l|a2  II  <  ■■■  <  Ikm-ill  <  \\am\\,  where  ajt  j  =  l,...,m 
are  columns  of  A.  Then  the  bounds  of  the  EGi(Q)  in  (70)  are 
tighter  than  those  in  (69). 


(\\am\\2 -\\amM\2)2' 

(80) 

To  compare  these  bounds  to  Rung  and  Suter  [9] ,  set  ||  c|| 2  = 
0,  and  the  bounds  for  EigGapj(Q)  in  (70)  become 

HGi(A)<EGi(Q)<HGi(A)HG2(A).  (81) 


+  ^(HG,(A)-1)HG2(A) 

A 


Under  these  conditions,  the  lower  bound  agrees  with  Rung 
and  Suter  while  the  upper  bound  is  tighter. 

Let  A  E  K"xm  be  an  ordered  hub  matrix.  Let  A  E 
Knxk  be  a  hub  matrix  obtained  by  removing  the  first  n  -  k 
nonhub  columns  of  A  with  the  smallest  Euclidean  norms. 
This  corresponds  to  the  MIMO  beamforming  problem  by 
comparing  k  of  m  transmitting  antennas  with  the  largest 
signal-to-noise  ratio  (see  [9]).  The  ratio  MQ)/M(Q)  with 
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Q  =  A‘A  describes  the  relative  performance  of  the  resulting 
systems.  It  was  shown  in  [9]  that  for  k  >  2 


h( 

112  "  WQ) 


I  dm  II  +11  am-\ II 


(82) 


II  Am  II  +11  1  I 

By  Theorems  1 1  and  12,  we  have 

f4(\\am-k+i\\2/\\am\\2)  <  h(Q)  U(\\cLm-A2/\\am\\2) 
/4(||flm-il|2/||flm||2)  “'UQ)-  /4(||fl1||2/||flm||2) 

(83) 

By  Lemma  18,  the  lower  and  upper  bounds  for  the  ratio 
Ait(Q)/Am(Q)  in  (83)  are  better  than  those  in  (82).  In 
particular,  when  k  =  1,  the  matrix  Q  corresponds  to  the 
hub,  as  such,  it  reduces  to  Q  =  [||am||2];  hence,  Ai(Q)  = 
||am||2.  Therefore,  an  estimate  of  the  quantity  ||flm||2/Am(Q) 
was  given  in  [9]  as  follows: 

ii  1 1 2  ii  1 1 2 


l|flml|2+l|flm-ll|2  A  m(Q) 
By  Theorems  1 1  and  12,  we  have 

II2  <  IlflJl2  < 


<  1. 


Clyy 


(84) 


(85) 


/4(lk„-ill2/lkJ2)  /4(llcnll2/llflmll2)' 

Again,  by  Lemma  18,  the  lower  and  upper  bounds  for  the 
ratio  l|flmll2/Am(Q)  in  (85)  are  better  than  those  in  (84).  We 
can  simply  view  (84)  and  (85)  as  degenerate  forms  of  (82) 
and  (83),  respectively. 

5.  Numerical  Examples 

In  this  section,  we  will  numerically  compare  the  lower 
and  upper  bounds  of  eigenvalues  for  arrowhead  matrices 
estimated  by  three  approaches.  The  first  approach  is  due  to 
Cauchy,  denoted  by  C  and  based  on  (2)-(4).  The  second 
approach,  denoted  by  SS,  is  based  on  eigenvalue  bounds 
provided  by  Theorems  11  and  12.  The  third  approach, 
denoted  by  WS,  is  based  on  Wolkowicz-Styan’s  lower  and 
upper  bounds  for  the  largest  and  smallest  eigenvalues  of  a 
symmetric  matrix  [15].  These  WS  bounds  are  given  by 


a  -  sp  <  Ai(Q)  <  a  — 
a  +  ^  <  A m(Q)  <  a  +  sp, 


(86) 


where  Q  e 


|mxm 
c2 


is  symmetric,  p  =  \!m-  1,  a  = 


trace(Q)/m,  ands2  =  trace(QfQ)/m  -  a2. 

Example  20.  Consider  the  directed  graph  in  Figure  1,  which 
might  be  used  to  represent  a  MIMO  communication  scheme. 

The  adjacency  matrix  for  the  directed  graph  is 


which  is  a  hub  matrix  with  the  right-most  column  cor¬ 
responding  to  node  4  as  the  hub  column.  The  associated 
arrowhead  matrix  Q  =  AT A  is 


Q  = 


10  0  1 
0  2  0  2 
0  0  11 
12  14 


(88) 


The  eigenvalues  of  Q  are  0,1,1.4384,5.5616.  Corollary  3 
implies  1  being  the  eigenvalue  of  Q.  By  Corollary  4,  the 
following  matrix 


Q  = 


1  0  V2 
0  2  2 
4  2  4 


(89) 


should  have  eigenvalues  of  Ai(Q)  =  0,  Xi{Q)  =  1.4384, 
Aj(Q)  =  5.5616.  The  C  bounds,  SS  bounds,  and  WS  bounds 
for  the  eigenvalues  of  the  matrix  Q  are  listed  in  Table  1 .  For 
Ai  (Q),  the  lower  SS  bound  is  best;  next  is  the  lower  C  bound 
followed  by  the  lower  WS  bound;  the  upper  SS  bound  is  best; 
next  is  the  upper  WS  bound  followed  by  the  upper  C  bound. 
For  A2(Q),  SS  bounds  and  C  bounds  are  the  same  as  the  C 
bounds.  For  A3(Q),  the  lower  SS  bound  is  best;  the  lower  C 
bound  and  WS  bound  are  the  same;  the  upper  SS  bound  is 
best;  next  is  the  upper  WS  bound  followed  by  the  upper  C 
bound.  In  conclusion,  the  SS  bounds  are  best. 

The  bounds  of  the  eigengap  of  Q  provided  by  (69)  are 


2  <  EGi  (Q)  <  6 
while  the  bounds  provided  by  (70)  are 

2.6861  <  EGi(Q)  <  5.6458. 


(90) 


(91) 


0  10  1' 

Example  21.  We  consider  an  arrowhea 

0  0  11 

‘  1  0  V2 

A  = 

,  (87) 

0  10  1 

Q  = 

0  6  2 

10  0  1 

i/2  2  7 

Therefore,  the  bounds  by  (70)  are  tighter  than  those  by  (69). 
This  numerically  justifies  Proposition  19. 


(92) 
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Table  1:  C  bounds,  SS  bounds,  and  WS  bounds  for  Example  20. 


Eigenvalues 

C  bounds 

SS  bounds 

WS  bounds 

Lower 

Upper 

Lower 

Upper 

Lower 

Upper 

Ai(Q)  =  0 

-0.4142 

1 

-0.3723 

0.3542 

-1 

0.6667 

A2(Q)  =  1.4384 

1 

2 

1 

2 

A3(Q)  =  5.5616 

4 

7.4142 

5.3723 

5.6458 

4 

5.6667 

Table  2:  C  bounds,  SS  bounds,  and  WS  bounds  for  Example  21. 

Eigenvalues 

C  bounds 

SS  bounds 

WS  bounds 

Lower 

Upper 

Lower 

Upper 

Lower 

Upper 

A  i(Q)  =  0.6435 

-0.4142 

1 

0.1270 

1 

0 

2.3333 

A2(Q)  =  4.6304 

1 

6 

4 

6 

A3(Q)  =  8.7261 

6 

10.4142 

7.8730 

9 

7 

9.3333 

Table  3:  C  bounds,  SS  bounds,  and  WS  bounds  for  Example  22. 


Eigenvalues 

Lower 

C  bounds 

Upper 

Lower 

SS  bounds 

Upper 

WS  bounds 

Lower 

Upper 

Ai(Q)  =  0.6192 

0.5000 

1 

0.2753 

0.8820 

0.2137 

0.9402 

A2(Q)  =  1.3183 

1 

2 

1 

2 

A3(Q)  =  3.0624 

2 

3.5000 

2.7247 

3.1180 

2.3931 

3.1196 

The  eigenvalues  of  Q  and  the  corresponding  C  bounds, 
SS  bounds,  and  WS  bounds  for  the  eigenvalues  of  Q  are  listed 
in  Table  2.  For  Ai  (Q),  the  lower  SS  bound  is  the  best;  next  is 
the  lower  WS  bound  followed  by  the  lower  C  bound;  the  SS 
upper  bound  and  the  C  upper  bound  are  the  same  and  are 
better  than  the  upper  WS  bound.  For  the  lower  SS 

bound  is  better  than  the  lower  C  bound;  the  upper  SS  bound 
is  the  same  as  the  upper  C  bound.  For  A3 (Q),  the  lower  SS 
bound  is  the  best;  next  is  the  lower  WS  bound  followed  by 
the  lower  C  bound;  the  upper  SS  bound  is  the  best;  next  is 
the  upper  WS  bound  followed  by  the  upper  C  bound. 

Example  22.  We  consider  an  arrowhead  matrix  Q  as  follows: 


Q 


0  2  1 


(93) 


6.  Conclusions 

Motivated  by  the  need  to  more  accurately  estimate  eigengaps 
of  the  system  matrices  associated  with  hub  matrices,  this 
paper  provides  an  efficient  way  to  estimate  the  lower  and 
upper  bounds  of  arrowhead  matrices.  Improved  lower  and 
upper  bounds  for  the  eigengaps  of  the  system  matrices  are 
developed.  We  applied  these  results  to  a  wireless  computation 
application,  and  subsequently  we  presented  several  numeri¬ 
cal  examples.  In  the  future,  we  will  plan  to  extend  our  results 
to  hub  dominant  matrices,  which  will  allow  hub  matrices 
with  correlated  nonhub  columns. 
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